Endpoint thermodynamics of an atomic Fermi gas subject to a Feshbach resonance 
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The entropy and kinetic, potential, and interaction energies of an atomic Fermi gas in a trap are 
studied under the assumption of thermal equilibrium for finite temperature. A Feshbach resonance 
can cause the fermions to pair into diatomic molecules. The entropy and energies of mixtures of 
such molecules with unpaired atoms are calculated, in relation to recent experiments on molecular 
Bose- Einstein condensates produced in this manner. It is shown that, starting with a Fermi gas 
of temperature T = O.IT^, where Tp is the non-interacting Fermi temperature, an extremely cold 
degenerate Fermi gas of temperature T < O.OlTp may be produced without further evaporative 
cooling. This requires adiabatic passage of the resonance, subsequent sudden removal of unpaired 
atoms, and adiabatic return. We also calculate the ratio of the interaction energy to the kinetic 
energy, a straightforward experimental signal which may be used to determine the temperature of 
the atoms and indicate condensation of the molecules. 

PACS numbers: 



I. INTRODUCTION 

The formation of cold diatomic molecules from an 
atomic Fermi gas offers one the opportunity to study the 
Bose-Einstein condensate to Bardeen-Cooper Schrieffer 
(BEC to BCS) crossover in an atomic system 
Cold atomic gases present several advantages in investi- 
gations of this long-standing problem: they are impurity- 
free; their temperature, density, and interaction strength 
can be tuned over wide ranges; and they are easily ma- 
nipulated with lasers. The key to obtaining supcrfluid 
states such as BEC and BCS is to cool the gas to a suf- 
ficiently low temperature Q. 

In recent experiments, Feshbach resonances have been 
used to transform a cold or deg enerate Fermi gas into 
a molecular BEC H 0, S H Efl 113 . A Feshbach reso- 
nance is a resonance between an unbound atomic state 
and a bound molecular state brought about by the Zee- 
man effect, in practice by changing a uniform magnetic 
field over the region of the gas [IJ ■ The molecular pair- 
ing mechanism is a few-body effect brought about by the 
Feshbach resonance. The BCS pairing mechanism is a 
many-body effect due to an instability at the Fermi sur- 
face. This interplay between few and many body physics 
is an outstanding problem. Various theories have been 
proposed to model the crossover region and molecule for- 
mation [IllllllllllllllllllllMillp. However, 
in order to study the crossover thoroughly it is desirable 
to have clearly defined and demonstrated endpoints. 

In the following, basic thermodynamic quantities re- 
lated to the endpoints of this problem are calculated: in 
particular, the kinetic energy, potential energy, interac- 
tion energy, and the entropy. Two cases are considered: 

(a) a degenerate atomic Fermi gas of two spin states in- 
teracting via a weakly attractive contact potential; and, 

(b) an atom-molecule mixture in which the atoms and 
molecules interact both with themselves and with each 
other via a repulsive contact potential, as was the case 
in the experiment of Ref. |^ . The atomic interaction pa- 



rameter is given by fc^|a|, where kp is the non- interacting 
Fermi momentum and a is the s-wave scattering length. 
The s-wave scattering length is a function of the magnetic 
field in the region of the resonance and typically changes 
over many orders of magnitude as well as in sign. Case 
(a) treats the parameter regime a < and kp\a\ ^ 1, for 
which a weakly interacting BCS-type transition is possi- 
ble. Case (b) treats the regime a > and kp\a\ <^ 1, in 
which case the molecules may Bose condense. In case (b) 
both a Boltzmann-like molecular gas of T/Tbec ^ 1 and 
a degenerate molecular gas of of T/Tbec *C 1 are consid- 
ered. Two applications which demonstrate the usefulness 
of these calculations in relation to present experiments 
are illustrated. 

Firstly, we present a straightforward signal for the on- 
set of condensation in the molecules. The ratio of the 
interaction energy to the kinetic energy is an established 
observable ■ The power law in this ratio changes as 
the molecules are condensed. Moreover, this quantity 
can be used to determine the temperature, which, par- 
ticularly in case (b) for T <^ Tbec, would be difficult to 
determine for the atoms from momentum distributions 
alone, since they arc pushed to the edges of the molecu- 
lar cloud by the molecular mean field, as we shall show. 
This ratio is thought to be a universal constant in the 
unitarity limit, i.e., when kp\a\ > 1, where it is called f3 
in the literature 24, 25] . 

Secondly, we propose a new way to obtain a sufficiently 
low temperature for a weakly interacting BCS transition. 
A recent proposal [2^ showed that adiabatic transition of 
the resonance from degenerate molecules to atoms leads 
to cooling. This follows from the fact that the entropy 
of an ideal fermionic gas S a. oc T while the entropy of an 
ideal bosonic gas Sra oc [23 . Holding the entropy con- 
stant results in a decreased temperature for sufficiently 
small T. The addition of weak interactions does not 
change this fact . In the experiment of Ref. 0| , it was 
found that, after tuning the atoms over the resonance so 
that they formed molecules, about 15% of the atoms re- 
mained unpaired [2^ . In such an atom-molecule mixture. 
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at low temperature, the majority of the entropy resides 
in the atoms. Therefore, suddenly removing the atoms is 
the same as removing entropy. The technique of remov- 
ing atoms while leaving the molecules in place has been 
experimentally demonstrated Adiabatically disso- 

ciating the molecules then leads to an extremely cold 
fcrmionic gas of temperatures T < O.OITf, where Tp is 
the Fermi temperature. We term this cooling method 
entropic cooling. Given the experimental limitations of 
evaporative cooling, which for a BEC is T ~ 0.25 Tbec 
and for a Fermi gas is T ~ 0.05 Tp (see, for example, 
Refs. llOl and 0,1^, respectively), entropic cooling pro- 
vides an attractive alternative. 

We note that a recent experiment has successfully tra- 
versed the crossover from a molecular BEC to an atomic 
Fermi gas in a sufficiently adiabatic manner to achieve 
entropic cooling |3^ . 

The presentation may be outlined as follows. In SccHTI 
the model is presented. In Sec. IIIII the thermodynamic 
quantities are calculated for a degenerate Fermi gas, in 
some detail. In Sec. IIVI the same quantitites are pre- 
sented for the atom-molecule mixture at temperatures 
both above and below Tbec- In Sec. [V] the two experi- 
mental applications are considered: a signal for the on- 
set of molecular BEC; and the use of entropic cooling to 
achieve temperatures sufficiently low for a BCS transi- 
tion. Finally, in Sec. I VII we conclude. 



II. THE MODEL 

Our model consists of the following assumptions. 
Firstly, it is supposed that for negative scattering length 
T > Tbcs and kp\a\ ^ 1. This assures that the system is 
far from resonance and that there is no superfluid phase. 
S econdly, for positive scattering length the dilute limit 
a/ nmafnm *C 1 is assumed for the molecules, where IS 
the average molecular density and the molecule-molecule 
and atom-molecule scattering lengths amm, Oam a, 
with a the s-wave atom-atom scattering length. The 
proportionality constants are taken to be 0.6 and 1.2, 
respectively, in accord with Petrov et al. |l6j| . where it 
is calculated explicitly in certain limits from the four 
body problem. Thirdly, it is assumed that the atoms 
and molecules can be treated as distinct entities: the 
former are fermions and the latter are bosons, whose 
thermodynamic properties can be calculated separately. 
Fourthly, it is assumed that the system is in thermal equi- 
librium. However, it is not necessary to require chemi- 
cal equilibrium. The former requires a sufhciently high 
rate of elastic two body collisions {A + A <-> A + A, 
M + A <-> M + A); as the latter is brought about by 
three-body collisions {A + A + A<-^M + A), it may be 
substantially slower. This indeed appears to be case in 
certain experiments 0, 0, |2^ [s^ . 

Fifthly, it is assumed that the atom-atom coupling 
for binary interactions is given by the s-wave limit g ~ 
2'!Th'^a/mr^ where = rna/'2 is the reduced mass. The 



atom-molecule and molecule-molecule coupling is then 
.9am = 0.9(7 and gmm = O.Sg according to Ref. [T6l |. 
Sixthly, it is assumed that both atoms and molecules 
are subject to the same isotropic harmonic trapping fre- 
quency uj. Thus their densities shall also be isotropic, 
as well as all other thermodynamic quantities [s^. Sev- 
enthly, it is assumed that for negative scattering length 
there are an equal number of spin up and spin down 
atoms. It is now possible to change the scattering length 
from negative to positive, selectively remove excess atoms 
of both spin states, and then change the scattering length 
back to negative, thereby resulting in the precise balance 
assumed here |29|. This results in a pleasant symme- 
try in the thermodynamic expressions, as well as avoid- 
ing Fulde-FerrcU-Larkin-Ovchinnikov phases when a BCS 
transition is sought (sec Ref. ^-nd references therein). 
Eighthly, the local density and semi-classical approxima- 
tions, which are valid for ksT 3> 1/p{Ef), are assumed 
to hold, where p{Ep) is the density of states at the Fermi 
surface. 

The critical temperature for condensation of an ideal 
gas shall be denoted as 
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where Nm is the number of molecules and C, is the Rie- 
mann Zeta function. The shift in this tem per ature caused 
by interactions may be calculated to be |3fil | 
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in a harmonic trap for the assumed regime of 
\/nmafnni < 1, whcrc 



(3) 



is the harmonic oscillator length. Therefore a ninth as- 
sumption is that this shift is small, so that it may be 
neglected in our calculations ^3 . This is consistent with 
the second assumption above. Typical shifts in experi- 
ments arc less than 10%. Similarly, the Fermi tempera- 
ture will be denoted by its ideal value of 
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where Na is the number of atoms. 

Lastly, for T <^ TbeCj the Thomas- Fermi profile shall 
be assumed for the bosonic molecules: 
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where the subscript m refers to molecules, to„j — 2ma 
and /i„j is the molecular chemical potential. Note that 



(r) = for r > i?„ 
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III. NEGATIVE SCATTERING LENGTH: 
ATOMS ONLY 

We first consider the case for which the scattering 
length is negative, so that there are no molecules. We 
calculate the entropy and the energy. The energy is di- 
vided into kinetic, potential, and interaction parts. For 
an isotropic harmonically trapped cold atomic gas, the 
potential energy 



V{r) = -TOaW^r^ 
represents the trap energy and 
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(7) 
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is the interaction or mean field energy, where na{r) is the 
total mean field particle density, i.e., including both spin 
states, and g is the interaction strength. The expressions 
are derived in some detail in this, the simplest case, in or- 
der to serve the reader as a model for similar calculations 
which are only sketched in the following sections. 



A. Energy 

The expressions for the energy in the local density ap- 
proximation arc given by integrals over phase space of 
the form 
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where 



v{r,p) = {cMP{E{r,p) - p)] + l]-^ 
is the mean occupation number, 
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and 



E{r,p) 



p. 



2mr 
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+ V{r) + C/,„f (r) 



(9) 
(10) 
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is the total energy in phase space. The factor of 2 in 
front of the integrals in Eqs. © and Hll|) is due to the 
number of spin states. 

One may obtain self-consistent expansions to first or- 
der in kp\a\ and second order in T/Tp, where 



(15) 



is the Fermi momentum defined with respect to an ideal 
gas. The key is to set na{r) = n^\r) in Unii{r) and 
Eq. l(TT)l . where 



(17) 



is the non-interacting density profile, which may be ob- 
tained by integrating Eq. (|12|l over momentum space. 
The function g„(x) in Eq. is the Bosc function [27l |. 
The details of the calculation are presented in App. El 
One finds 
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These expressions may be shown to be self-consistent 
in two ways. Firstly, one observes that the virial theorem 



2£'kin — 2i?, 



■pot 



3Si„t = 



(21) 



holds up to first order in k^\a\ and to second order in 
T/Tp. Secondly, note that the total energy at zero tem- 
perature for an ideal gas is i^^io' = -£^kin+-Epot = ^NEp, 
as may also be calculated directly from the density of 
states for an ideal Fermi gas in a harmonic trap. The ex- 
pansions given by Eqs. lfTH|) . ((T^ . and are accurate 
to 1%, 5%, and 24% for k°p\a\ = 0.1 and T/T° = 0.1. 
For higher temperature or stronger interaction strength, 
the accuracy becomes quite poor, for example 72%, 35%, 
and 233% for k°p\a\ = 0.5 and T/T^ = 0.5. In this 
case it is preferable to evaluate the energy integrals self- 
consistently in the mean field by recursive use of Eqs. ©- 
(|ll|l with the density profile 



n,(r) = 2 



d^p 
{2Trhy 



■v{r,p) 



(22) 



while holding the total number of atoms, given by the 
integration of Eq. H22|l over volume, constant. This gives 
the exact result, which we have used to evaluate the ac- 
curacy of the expansions. 

In Fig. n is shown the ratio of the interaction en- 
ergy to the kinetic energy in the degenerate regime. 
This quantity is a straightforward experimental observ- 
able |23j which can be used to determine the tempera- 
ture of the system. The figure shows the self-consistent 
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FIG. 1: Shown is the ratio of the interaction energy to the 
kinetic energy as a function of T/Tp for a weakly interacting 
degenerate Fer mi g as. This ratio is an established experimen- 
tal observable [2al that may be used to measure the degen- 
eracy of the gas. This ratio, called /3 in the unitarity limit, 
i.e., when kp\a\_>l, is thought to be a universal constant 
in that regime |24l . I2R| | . Here, the temperature dependence of 
the ratio is illustrated for increasing values of fc^|a|. 



For an ideal Fermi gas in a harmonic potential, the in- 
tegral of Eq. (|25|l can be performed exactly (see Eq. (|A14(I 
in Add. One obtains 
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In the degenerate limit, f3^a S> 1, the Bose functions may 
be expanded to obtain (see Eq. IjAlip in App. ^ 



S = ksn^N 



T 
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to lowest order in T/T^. As in Sec. 111111 corrective 
terms due to the interactions may be obtained in powers 
of fcj!n|a|. The method for doing so is outlined in Add. 1X1 
The result is 



(28) 



to lowest order in T/Tp. For values of k'p\a\ < 1/2, as 
we have assumed, the correction is less than 10%. 



solution which makes no assumptions on the smallness 
of the mean field or the temperature. This is partic- 
ularly important for the upper curve, which treats the 
almost strongly interacting regime. In a recent work, it 
was shown that values of fc^|a| likely to lead to a BCS 
transition are on the order of one- half ^9 . 



B. Entropy 

The general combinatorial expression for entropy is 

(23) 



S = -ks^px Inpx 

{A} 



where {A} represents all states of the system and p\ is 
the probability to be in a given state [s^ . For fermions in 
the local density approximation in an isotropic system, 

S = -2kB^[i^\niy+{l-,^)\n{l-iy)], (24) 

where the factor of 2 is due to the two atomic spin states 
and ly is defined in Eq. (|12|l . Equation (|24|l can be inte- 
grated over momentum to obtain 
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The details are given in App. This is the general 
expression for the entropy of the fermions. 



IV. POSITIVE SCATTERING LENGTH: 
MOLECULE-ATOM MIX 

We next consider the case of positive scattering length, 
for which there is a mixture of atoms and molecules, as 
in the experiment of Ref. In the following, it shall 
be assumed that the molecular mean field acts on the 
atoms, but that the effect of the atomic mean field on 
the molecules is negligible. This is justified for small 
fractions of atoms, which is the case experimentally. The 
effect of the atomic and molecular mean fields on them- 
selves, respectively, shall be calculated to first order in 
the interactions, as was done in the previous section for 
negative scattering length. 



A. Above the Molecular Condensation 
Temperature 

1. Molecules 

The entropy and kinetic, potential, and interaction 
energies of the molecules may be calculated by similar 
methods to those presented explicitly in Sec. lIIII The lo- 
cal density approximation for the energy in phase space 
is 



E{r,p) 
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One may take the mean field density profile to be the non- 



interacting one, i.e., n^, 



1^(0) 



r(r) (as in Eq. ^) 
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with all a subscripts changed to m, where the factor of 
one half is to account for the number of spin states. Then 
the resulting energies arc: 
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to cubic order in Tb^q/T and leading order in the mean 
field in Eq. where 
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FIG. 2: Shown is the density profile of a Fermi gas (solid 
curve) in the presence of a molecular condensate (dashed 
curve), the latter being in the Thomas- Fermi limit with 
the Thomas- Fermi radius. The atoms occupy a spherical shell 
outside the molecular region. Note that the parameters of the 
JILA experiment have been used for T/TgEc = 0.25 with 
the molecule-molecule and atom-molecule scattering lengths 
derived in Ref . . 



is the approximate binding energy of the molecule. Ap- 
pendixElgives the perturbative expression for the molec- 
ular fugacity, which is useful in obtaining the energies, 
in that one must eliminate the chemical potential in or- 
der to obtain energies depending on temperature and 
atomic/molecular number alone. 

Note that the large temperature non-interacting limit 



is given by E^ 
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kin 



■E, 



pot 



SN^kBT. This agrees 



with a simple calculation based on the density of states: 
E^^t = SdEE^eP^>'-'-^y{nLuf = 3iV„,fcBT + 0(T3). 
The expansions given by Eqs. H3U|I . (|31|l . and (|32|l are 



accurate to 0.66%, 2.4%, and 6.4% for fc? 



0.5 and 



T/TggQ ~ 1. For higher temperature or weaker interac- 
tion strength, they are accurate to better than 1%. 
The physically measurable quantity 
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(35) 



to leading order, in contrast to the power law for Bose 
condensed molecules, as shaU be discussed in Sec. IIVBI 
The entropy is calculable by similar techniques as were 
used to obtain the energies, either from the combinatoric 
expression similar to the calculation of Sec. lIIIBl or from 
the grand potential (see, for example, Ref. |2a|). 



Note the inclusion of the molecular mean field. The 
atom-atom and atom-molecule interaction energies are 
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respectively, while the kinetic and potential energies are 
calculated in analogy with Eqs. Q and (|10|l . One finds 
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2. Atoms 



The energy of the atoms in phase space is 
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where the expansions have been made to cubic order 
in the small parameter T^^^/T and to lowest order 
in the atomic and molecular mean fields in Eq. H36|l . 
Appendix [0 gives the perturbative expressions for the 
atomic and molecular fugacities, which are useful in 
obtaining the energies. Note that, as in Sec. IIV A II 
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the large temperature non-interacting limit is given by 
-E'tot = £'kin + ^'pot = SNakgT. The expansions given 
by Eqs. ij^ll, (EU and ^ are accurate to 0.006%, 
1%, 4%, and 4% for fc^|a| = 0.5 and T/T^ec = 2- For 
higher temperature or weaker interaction strength, the 
accuracy rapidly improves. In the most extreme limits, 
e.g. I^/TggQ = 1, the accuracies arc on the order of 10% 
to 50%. 

The correct virial theorem must be modified from that 
of Eq. (EB: 



2£'kii 



2E, 



pot 



o iT'aa 



(43) 



Note that, according to the above expressions, the con- 
dition for chemical equilibrium, 2fj,a = —Eb + fJ-rm which 
we have not assumed here jT3|, would give 



/rpO \ 3/2 



(44) 



for T^^JT « 1. 

The entropy of the atoms may be calculated by inte- 
grating numerically over the general expression for the 
entropy density given by Eq. (|25|l with 



B. Below the Molecular Condensation 
Temperature 

1. Molecules 

The chemical potential of the molecules may be calcu- 
lated from the Thomas- Fermi profile rv^f{r) to be 



3-5 



2/5 



(.9mmW™)'/'(™rnCc>2)3/^ (46) 



The energy in phase space is given by the energy of 
Boguliubov quasiparticle excitations (see and refer- 
ences therein): 
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The energies may then be solved for as 
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An expansion may also be developed in order to deter- 
mine a perturbative expression in the Boltzmann limit in 
the same fashion as was done for the energies. 
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where the integrals for the kinetic, potential, and inter- 
action energies are defined respectively by 
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Note that the virial theorem, Eq. H21|l . holds. 

The entropy has been calculated elsewhere in detail 
from the exact density of states |2^ : 

= [filJiTTn^oj^)]{2V2zt&nh~^[V2l/{l + z)] 
+4:z^/^ - V2z^[n + 2 tan"^((l - z)/V2i)] 

+ (l + z)2[0o-sin(40o)/4]}, (52) 

where z = e/ is the rescalcd energy and ^o = 
cos~-'-(l/Vl + z). The result is 
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FIG. 3: Atomic Fermi gas. Shown is the ratio of the in- 
teraction energy to the kinetic energy as a hinction of the 
degeneracy T/Tp, all in the presence of the molecular mean 
field. Since the atoms are pushed to the outside of the 
molecular mean field in space, the degeneracy cannot be ex- 
tracted in a simple way from the momentum distribution; 
this plot serves as an alternative method to determine their 
degeneracy. The left hand curve is calculated for a Bose- 
condensed molecular mean field, i.e., T ^ T^ecj while the 
right hand one is calculated in the opposite limit. Note that 
r]J = 1.53(iVa/iVm)^''^TiEc; for the parameters shown here, 

= 1.11 r^Ec- 



2. Atoms 

The energy of the atoms in phase space is 
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Inserting the Thomas-Fermi profile for the molecular den- 
sity, one finds 
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(57) 

Thus the atoms experience an expulsive harmonic poten- 
tial |39j due to the molecular mean field and are pushed 
to the outside of the molecular cloud. The effect is illus- 
trated in Fig. 121 To calculate the atomic density profile 
to first order in the interactions, one follows similar pro- 
cedures as outlined in previous sections. The result for 
(r) is given by Eq. with 



/3 {^ia + Imoj'^r^ - S^i™) |r| < i?„ 



/3 (/ia — hmuj'^r'^^ 



\r\ > Rr. 



Unlike in the other cases studied for atoms, we did 
not find a simple way to expand the Bose functions in 



FIG. 4: Molecular Bose gas. Shown is the ratio of the inter- 
action energy to the kinetic energy as a function of T/T^ec- 
This ratio may be used as a signal of condensation since, 
for the molecules, the power law of £int/Ski„ oc {T/TiEcY 
changes for 5 ~ —5/2 above the phase transition to J ~ — 1 
below. The left hand curve is calculated for a Bose-condensed 
molecular mean field, i.e., T ^ Tbeci while the right hand 
one is calculated in the opposite limit. The parameters of the 
JILA experiment have been used for plotting purposes 0. 



the integrands of the energies. Therefore the energies 
are left determined up to an integration which may be 
performed numerically for a given parameter set, as in 
Eqs. 63 - (EB. 

The entropy of the atoms may again be calculated by 
integrating numerically over the general expression for 
the entropy density given by Eq. H25|) with 



1 



2 2 



1 



gn-air) 



TF 



(r)] (59) 



V. APPLICATION TO EXPERIMENTS 

The entropies and energies calculated in Sees. IIIII 
and IIVI are useful in understanding experiments. In the 
following, two examples are given. In Figs. I^lsl and Ta- 
ble m the parameters of the recent JILA molecular BEC 
experiment |^ are used throughout. 



A. Evidence of molecular condensation and 
atom-molecule temperature 

The dependence of the ratio of the interaction energy 
to the kinetic energy on temperature is a useful experi- 
mental observable. It is obtained via time-of-flight mea- 
surements with and without interactions; the interactions 
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can be rapidly switched on or off via a Feshbach reso- 
nance jl^. This ratio allows one to determine the de- 
generacy in the case of negative scattering length (see 
Fig. Q). In the case of positive scattering length, when 
one has an atom-molecule mixture, it can again be used 
to calculate the temperature of the atoms. This is partic- 
ularly important as their temperature dependence can- 
not be determined in a simple way from their momentum 
distribution, due to the fact that the atoms occupy a 
spherical shell outside the region of the molecular mean 
field (see Fig. In Fig. El is shown this ratio, calcu- 
lated for T both above and below TggQ, from the inte- 
gral equations of Sees. IIV A 21 and HV B 21 In the case of 
molecules, the dependence of the ratio of interaction to 
kinetic energy on degeneracy signals the onset of conden- 
sation. Above T^^^^ E,^t/E^^ cx {T /T^^cV'/'^ while 
below i?int/£^kin fx {T /T^-qq)^^ . The result is illustrated 
in Fig. 21 as calculated from the energy equations of 
Sees. IIV A II and IIV B II Both figures use the parame- 
ters of the JILA experiment |^. Note that the fact that 
the left hand and right hand curves in Figs. and ^ do 
not meet is due to our use of the Thomas-Fermi approx- 
imation for the mean field of the molecular condensate. 
The use of this approximation changes the form of the 
left hand curves as T — > Tg^^ from below; the right hand 
curves are exact. 



B. A New Entropic Cooling Method 

The BCS transition temperature is given by 

Tbcs/T]^ = 0.277 exp(-^). (60) 

For k^p\a\ = 1/2, the temperature necessary to achieve 
BCS is Tbcs = 0.012 r]J. Present evaporative coohng 
methods have been unable to achieve temperatures below 
about T = O.OSTf 0|. Recently, it was suggested that 
adiabatic transition across a Feshbach resonance could be 
used to transform a molecular BEC into a very cold Fermi 
gas, of sufficiently temperatures to achieve BCS j26j . 
This is due to the fact that the entropy of fermions is 
proportional to T while that of bosons is proportional 
to T^. This scheme requires evaporative cooling of the 
molecular BEC. The typical maximal degeneracy achiev- 
able in an atomic BEC experiment is T/T^^^^ ~ 1/4. 

Consideration of an atom-molecule mixture presents 
another alternative. Rather than trying to lower the 
temperature of the system by evaporative cooling, one 
may instead attempt to decrease the entropy. Since en- 
tropy is held constant across the Feshbach resonance, this 
amounts to decreasing the final temperature of the Fermi 
gas. In Fig. \S[ is shown the total entropy of an atom- 
molecule mixture for the JILA parameters |^ according 
to the expressions derived in the preceeding sections. As 
illustrated in the figure, the atoms contribute the major- 
ity of the total entropy when the Bose gas is degener- 
ate. It has been demonstrated experimentally that the 




FIG. 5: Shown is the total entropy of an atom-molecule mix- 
ture (solid curve), as well as the separate contributions of the 
atoms (plus signs) and molecules (open circles). Clearly, for 
low temperature, the atoms carry the majority of the entropy, 
which is what makes the entropic cooling scheme outlined in 
Table effective. Note that the parameters of the JILA ex- 
periment have been used. 



atoms can be selectively removed from an atom-molecule 
mixture . We therefore propose the following cooling 
scheme: 

1. A pseudo-spin- 1/2 Fermi gas is evaporatively 
cooled to degenerate temperatures for a small neg- 
ative scattering length (fc^|a| ^ 1). 

2. A uniform external magnetic field is changed across 
a Feshbach resonance so as to cause the atoms to 
be transformed into molecules, slow enough to en- 
sure thermal equilibrium. The scattering length is 
now small and positive (fc^|a| ^ 1, or, equivalently, 
n^/^a < 1). 

3. A small number of atoms remain unpaired 0|. 
These atoms are suddenly removed from the sys- 
tem, thereby reducing the entropy. 

4. The system is swept adiabatically back across the 
Feshbach resonance to a small, negative scattering 
length. All molecules dissociate into atoms. The 
resulting Fermi gas has a reduced temperature. 

5. The sequence is repeated until temperatures needed 
to obtain a BCS transition are achieved. 

In TablcQlis presented a realization of this scheme, utiliz- 
ing parameters similar to those of Ref. @ , where typically 
85% of the atoms are converted into molecules. After 
three sweeps through the resonance, the final tempera- 
ture of T = O.OOlOTp is so low as to be out of the range of 
the semi-classical approximation. Detailed requirements 
on the adiabatic transition times and loss rates in the 
system have been presented elsewhere |26| . 
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TABLE I: An application of entropic cooling. The left two 
columns depict the degeneracy and total number of fermionic 
atoms for negative scattering length. After adiabatic tuning 
via a Feshbach resonance to positive scattering length (right 
facing arrow) , 84 % of the atoms are transformed into bosonic 
molecules [3. In the two right hand columns are shown the 
resulting entropy and degeneracy. The remaining atoms are 
then suddenly removed 1291 , thereby decreasing the entropy of 
the system, as shown in the third column of the next row. Af- 
ter tuning back (left hand arrow), the temperature is reduced. 
After several repetitions of this process, the degeneracy is high 
enough to achieve a weakly interacting BCS transition. 



T/Tf 


Na 


Switching 


Entropy (/cb) 


T/Tbec 


0.100 


200,000 




217,000 


0.640 


0.0504 


167,000 




91,300 


0.640 


0.0504 


167,000 




91,300 


0.380 


0.0147 


139,000 




22,300 


0.380 


0.0147 


139,000 




22,300 


0.120 


0.00103 


116,000 




1,301 


0.120 



corresponds to a large reduction in entropy, which, af- 
ter return to the negative scattering length, results in a 
decreased temperature. For example, for a conversion 
efficiency of atoms to molecules of 85%, as in Ref. @, 
and starting with a Fermi gas of degeneracy T /Tp = 0.1, 
after three sweeps the degeneracy is increased by two or- 
ders of magnitude to TjTp = 0.001 with a loss of less 
than half of the atoms plj . 

We note that entropic cooling was used successfully for 
the first time in a recent experiment |32l| . 

We thank Deborah Jin, Marcus Greiner, Cindy Re- 
gal, Gora Shlyapnikov, and Yvan Castin, for useful dis- 
cussions. We are grateful to Jochen Wachter, Jinx 
Cooper, and Simon Gardiner for a careful reading of 
the manuscript. The support of the U.S. Depart- 
ment of Energy, OfRcc of Basic Energy Sciences via 
the Chemical Sciences, Geosciences and Biosciences Di- 
vision is acknowledged. L. D. Carr also acknowledges 
the support of the National Science Foundation via grant 
no. jmps-drf 0104447| 



We note that, for the purposes of the table, it was 
assumed that ^am = 0.95 ^-nd g-mm = O.Sg, where g^m and 
5mm ai'c the couplings for atom-molecule and molecule- 
molecule interactions [T^ . This was necessary in order 
to properly treat the effects of the atomic and molecular 
mean fields. 



VI. CONCLUSIONS 

The entropy and kinetic, interaction, and potential en- 
ergy of fermionic atoms and bosonic molecules were cal- 
culated for temperatures both above and below Tp and 
^BEC ■ effect of the atomic and molecular mean fields 
were taken into account both perturbatively in analytic 
expansions and sclf-consistently, as illustrated in the fig- 
ures. Below the condensation temperature, the molecu- 
lar mean field was assumed to be in the Thomas-Fermi 
limit and the thermal cloud was neglected. The general 
method of performing these thermodynamic calculations 
was outlined. 

Two applications of the calculations were then pro- 
posed. Firstly, it was shown that the ratio of the kinetic 
to interaction energy is sufficient to determine the degen- 
eracy of the atoms. The same ratio, for the molecules, 
may be used to indicate the onset of Bose condensation, 
as the power law dependence on the degeneracy changes 
from an exponent of —5/2 to —1. Secondly, it was sug- 
gested that an adaption of entropic cooling [2^ could 
be used to cool a Fermi gas to sufficiently high degener- 
acy so as to achieve a weakly interacting BCS transition, 
without the need for evaporative cooling of either the 
fermionic atoms or the bosonic molecules. This scheme 
involves tuning adiabatically back and forth through a 
Feshbach resonance: on the positive scattering length 
side, the sudden removal of remaining unpaired atoms 



APPENDIX A: CORRECTION TO THE 
KINETIC, POTENTIAL, AND INTERACTION 
ENERGIES FOR A WEAKLY INTERACTING 
FERMI GAS 

Beginning with the general expressions for the energies 
given by Eqs. Q - (|ll|l . one may integrate over momen- 
tum to obtain 



Jo 

,3/2, 



^pot = Apot / dy2/'^'53/2(-e") 





Eint = ^int / dyy^'"^ 
Jo 



where 



u = (3 



and 



18N f T 



Akin - -^[rfo\ ^bT, 

12N ( T\ 
Apot — — I I , 



6N / T 



7/2 



(Al) 
(A2) 

(A3) 

(A4) 
(A5) 

(A6) 
(A7) 



kBTk"p\a\. (A8) 



To resolve the integrals to first order in the interac- 
tions, it is sufficient to assume a mean field given by 
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the non-interacting density profile of Eq. (|16() . Then, 
one may expand the Bose functions in the integrands of 
Eqs. (EU) - as 



gn {- exp [f3na 



y - "53/2 



')]} 



- gn[~ exp(/3^a - y)] 
-"53/2 [- exp(/3^a - 2/)]5ri-i [" exp(/3^a - ?/)] , (A9) 



where 



2m 



(AID) 



is the small expansion parameter. The maximum value of 
the Bose function in the non-interacting density occurs 
for y = for positive I3fj,a- In the limit /3/Za ^ 1, the 
Bose function may be expanded as [s^ 



9n i~en = - 



Tin) 



1 tt\ ,1 
- + — n-1 — 
no 



-I^(n-l)(n-2)(n-3)^ 
360 



Thus, in the degenerate regime, 



[ff3/2 



)] 



3/2 



(All) 



(A12) 



One may therefore take as a condition for the expansion 
of Eq. (ttflll : 



(A13) 



where Eb = h'^/{maa?) is the approximate binding en- 
ergy of the molecule. 

Using the integral relations 



9i(z) = 
93iz) 



2 

JO 

2 



dyy^^^g5/2ize ^) , 
dyy^^^g3/2{ze~y) , 



(A14) 



eic, one may resolve the first term in the expansion of 
Eq. ljA9p exactly. An additional useful relation, neces- 
sary for calculating the potential energy by integration 
by parts, is 



dygn {qe ^) = -.g„+i {qe . 



(A15) 



To resolve the order a term of Eq. (jASJ, one notes 
that there are two regimes of integration: (a) y <JC f3fia 
and (b) y ^ Pfia- For a degenerate Fermi gas, so that 
f3^a 3> 1, the two leading contributions to the integrals 
are obtained in regime (a). Therefore, in order to find the 
result to order (T/T^)^, for which the two leading order 
terms are required, it suffices to perform the integral by 
the use of the expansion given in Eq. (|Alip . Expanding 



around small y and integrating, one obtains a final result 
similar to Eqs. (|18|) - (|20|l but dependent on the chemical 
potential Ha- 

To obtain a self-consistent correction, the relation be- 
tween the chemical potential and the Fermi temperature 
can be calculated from the equation for the total number 
of fermions 



N = n, 



tfirdPp 



(A16) 



where rig = 2 is the number of spin states. One finds, by 
the same methods as outlined above, that 



= 1- 



^ 16 , n I 

\ fcp a 

3 315 ^' 



^ ) (A17) 



to second order in T/Tp and first order in /c^jaj. Elim- 
inating the chemical potential from the equations for 
the energies via substitution of Eq. (|A17|) . one obtains 
Eqs. (O - (HOI). 



APPENDIX B: CORRECTION TO THE 
ENTROPY FOR A WEAKLY INTERACTING 
FERMI GAS 

Equation H24|l may be conveniently rearranged as 

{r,p} 



where 



(B2) 



Assuming the semiclassical local density approximation, 
one can make the change of variables 

u = P[fia-Vir)~Urajir)] , (B3) 

„2 



p 



2ma 



S = j d^rS{r) 



(B4) 
(B5) 



Then, integrating the second term in Eq. ljBl|l by parts, 
one finds 

Sir) = A / dzz^/^- — (B6) 



A 



ks f 2™a 

2^ \j3r? 



1 + exp(2: — u) 

3/2 



(B7) 



This integral may be resolved as a sum of two Bose func- 
tions |23| of form gn{x), as given in Eq. H25|) . 
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APPENDIX C: DEPENDENCE OF THE 
FUGACITY ON TEMPERATURE IN THE 
BOLTZMANN LIMIT 



In order to calculate the thermodynamic quantities for 
the atom-molecule mixture at T > T^^^ in Sees. IIV A II 
and II V A2l the following pcrturbative expansions of the 
fugacities are useful. Starting with Eq. ljA16|l . and in- 
cluding the mean fields, the relations are 



C(3) 



f '^BEC ^ 




f ^BEC \ 









11/2 



[C(3)]^ 



^BEC 
T 



(CI) 



for the molecular fugacity and 




for the atomic fugacity, where 
Eq. 



11/2 



(C2) 



7 oc kp\a\ is defined in 
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